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Abstract 

The long-range quadrupole-quadrupole (~ R^^) and leading dispersion (~ R^^) interactions be- 
tween all pairs of excited Hg(6s6p) ^Pq, ^Pi, ^P2i s^^nd ^Pi atoms are determined. The quadrupole 
moments are calculated using the ab initio relativistic configuration-interaction method coupled 
with many-body perturbation theory. The van der Waals coefficients are approximated using pre- 
viously calculated static polarizabilities and expressions for the dispersion energy that are validated 
with similar systems. The long-range interactions are critical for associative ionization in thermal 
and cold collisions, and are found to be quite different for different pairs of interacting states. Based 
on this knowledge and the short-range parts of previously calculated potential curves, improved 
estimates of the chemi-ionization cross sections are obtained. 
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I. INTRODUCTION 



Potential-energy curves for the interaction of two mercury atoms in their lowest ex- 
cited configuration (6s6p) have previously been calculated and used to estimate the chemi- 
ionization cross sections for collisions of two such atoms However, those potential curves 
were not accurate at distances beyond the valence region, which are of great importance in 
thermal and cold collisions. Reactions at such energies are pertinent to modeling of fluores- 
cent lamps as well as recent considerations of mercury as a candidate for an atomic clock 
[2I. In the present work, the long-range interactions are calculated and combined with the 
short-rauge parts of the previously ealeulated poterrfal eurves H „r order to provide better 
estimates of the chemi-ionization cross sections. 

Chemi-ionization is allowed if the potential energy of the two interacting Hg atoms 
lies higher than the potential energy of ground-state Hg^ at any internuclear distance 
R. For Hg(3Po)+Hg(3Po), Hg(3Pi)+Hg(3Po), Hg(3Pi)+Hg(3Pi), Hg(3p2)+Hg(3Po), and 
Hg(^P2)-fHg(^Pi) only associative ionization 

Hg* (6s6p) + Hg* (6s6p) ^ Hg+ + , (1) 

utilizing the attractive interaction energy between Hg"*" and ground-state Hg, is possible at 
very low energies since their total excitation energies are smaller than the ionization potential 
of Hg. However, as discussed in Ref. based on the potential curves of Hgg and Hg2 , 
essentially all chemi-ionization in mercury at very low collision energies, even for the other 
reactant pairs, Hg(3p2)+Hg(3p2), Hg(iPi)+Hg(3Po), Hg(iPi)+Hg(3Pi), Hg(3p2)+Hg(iPi), 
and Hg(^Pi)+Hg(^Pi), can be expected to be of the associative type and will be referred to 
as associative ionization (Al) henceforth. 

For nonsymmetric atomic states (i.e., both atoms having J 7^ 0), the leading term in the 
long-range interaction is generally of the form C^/R^, where C5 can be positive or negative. 
However, we will show that J 7^ is neither a necessary nor sufficient condition for C5 7^ 0. 
If C5 is positive and sufficiently large, then reactions such as AI, which occur at short range, 
are precluded. In all cases there is also an attractive van der Waals interaction of the form 
Ce/R'. 

For the ten possible combinations of the ^Pq, ^Pi, ^P2, and ^Pi states that may re- 
act, there are total of 90 potential-energy curves describing the interactions. The AI cross 
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section of only one of these pairs — Hg(^Pi)+Hg(^Po) — has been measured experimentally 
(although one experiment apparently misidentified these reactants as Hg(^Po)+Hg(^Po))- 
This experimental value is particularly useful in validating the present cross sections since 
the autoionization width has not yet been calculated but is assumed to be large. The ex- 
perimental cross section, combined with the theoretical potential curves, suggests that the 
autoionization width is indeed sufficiently large to saturate ionization when the continuum 
is penetrated at short range. 

II. QUADRUPOLE MOMENTS OF ATOMS AND MOLECULAR C5 COEFFI- 
CIENTS 

The quadrupole moments of the mercury atom in its 6sQp ^Pi, ^P2, and ^Pi states (the 
quadrupole moment of the ^Pq atom vanishes by symmetry), as well as the off-diagonal ma- 
trix elements, are calculated using the ab initio relativistic valence configuration-interaction 
(CI) method based on the Brueckner orbitals [sl, which effectively includes core polarization. 
In addition, the random-phase-approximation (RPA) chain of diagrams is included. This 
method expands the two-particle basis functions as 

^(vr, J,M) = J2 Ck,i^k,ii7r, J, M), (2) 

k>l 

where vr is the parity of the state \1/ and J and M are the total angular momentum and its 
projection. The weights Ck,i are determined by solving the eigenvalue problem based on the 
Hamiltonian in the model space spanned by the basis functions defined in the subspace of 
virtual orbitals 

$fe/7r, J, M) = rik,l J2 4X.;..,-,«{n„.„m,}«{n„K„,n,}|0core) • (3) 

Here the sets {n, k, m} enumerate quantum numbers, the rj^^i are normalization factors, the 
are creation operators, and the quasivacuum state |Ocore) corresponds to a closed-shell 
core containing 78 of the 80 electrons of the Hg atom. The one-particle orbitals (pk are 
determined by solving a Dirac equation 

{ho + ^DHF + T.)(t>k = ek4>k ■ (4) 

Here ho includes the rest mass term, kinetic energy, and the Coulomb interaction with 
the nucleus. I^dhf is the Dirac-Hartree-Fock (DHF) potential due to core electrons. E is 
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the self-energy operator computed in the second-order approximation (see, e.g., Ref. {4]). 
Quahtatively, this correction describes a response of a valence electron to polarization of the 
core by the electron's own field and is the dominant correlation effect. The resulting one- 
particle orbitals are usually referred to as Brueckner orbitals. Once the Brueckner orbitals 
are computed, we solve the CI problem in the model space. The model space Hamiltonian 
includes the Coulomb interaction between valence electrons. We do not include so-called 



screening corrections to the Hamiltonian. As demonstrated {4], these corrections, although 
computationally expensive, are relatively insignificant for divalent atoms. 

With the determined CI wave functions we form a matrix element of the quadrupole 
operator. The associated one-particle reduced quadrupole matrix elements are given by 

/■oo 

= ('^.llC^'^ll'^i) / r2[G,(r)G,(r) + F,(r)F,(r)]dr, (5) 
Jo 

where C*^^-* is the normalized spherical harmonic [sl, and G and F are the large and small 
radial components of the relativistic wave function. Further we employ the RPA approxi- 
mation 6|], which accounts for shielding of an externally applied field by the core electrons 
and substitutes dressed matrix elements for the bare one-particle matrix elements. 

The quadrupole-quadrupole interaction energy between two atoms is given by [3, S Is] 

1 2 

^W = -^ (6) 

where 

= (2-,)^(2 + ,)r « 

The quadrupole spherical tensor is defined by^ 

i 

where the sum goes over all atomic electrons (though in the case of IIg(6s6p) only the 
unpaired 6p electron contributes). The quadrupole moments Q(^^'^^Pj) are defined, as 
conventional, by the component in the "stretched" M = J state: 

QC^'Pj) = 2C'^'Pj,j \Qo\ ''^'Pj,j) (9) 



^ Atomic units are used except where explicitly indicated otherwise. 
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(note the factor of 2; some definitions don't include this factor [10]) and is related to the 
reduced matrix element, given by Eq. ([5]), by a factor derived from the Wigner-Eckart 
theorem 



. ^ 2 J , 

cwe(J) = , 10 

^ -J J 

so that 

QC'^'Pj) = 2cwE(J)r+^Pj IIQII ''^'Pj). (11) 

Both the reduced matrix elements and the conventional quadrupole moments are given in 
Table [H The accuracy of these moments is expected to be a few percent (for example, we find 
that the RPA shielding reduces the CI values for quadrupole moments by less than 3%). We 
also note that nonrelativistically the ratio P2) / Pi) = —2, while in our calculations 
this ratio is —2.45. This deviation reflects the importance of relativistic corrections for the 
heavy mercury atom. 

For a particular pair of atomic states ^'^"^^Pj^ and '^^'''^^Pj^, the Hamiltonian H for the 
reduced quadrupole moment can be obtained in a double atomic basis. For two atoms in 
the same state, ^ = \Sa, La, Ja, Ma)\Sa, La, Ja, Mb) and V'' = \Sa, La, Ja, M'^)\Sa, La, Ja, M[) 
for the initial and final state, respectively. For the two like atoms in different states, 
{Sa, La, Ja) 7^ {Sh, Lf,, Jb), Symmetric and antisymmetric linear combinations of the func- 
tions on centers I and II must be formed, 

7/> = 2-'/\\Sa, La, Ja, Ma)l\Sb, Lb, Jb, Mb) H ± \Sa, La, Ja, Ma)n\Sb, Lb, Jb, Mb) l) (12) 

and 

ij' = 2-'/\\Sa, La, Ja, M'^)j\Sb, Lb, Jb, M',) n ± \Sa, La, Ja, M'^) n\Sb, Lb, Jb, Ml)i). (13) 

The quadrupole-quadrupole interaction is then given by 

{^\Q,Q.,W) = 

$:^(/.)(-i)^»--"(-i)^--^ f J. 2 j.w J. 2 Jb\ ^^iiQii^^^^iiQii,^ 

Y \'Ma Ml^J \-Mb -/i M^J 

(14) 
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in terms of the Wigner 3-j symbols and the reduced matrix elements, which satisfy 
(a||(5||&) = {—l)'^''~'^''{b\\Q\\a). By the properties of the 3-j symbols, we can get two se- 
lection rules: (i) the first (diagonal) term vanishes if < 2 or if J;, < 2 and (ii) the second 
(off-diagonal) term vanishes if \ Ja + Jb\ < 2. Thus the Hg('^Po)+Hg('^P2) collision possesses 
a long-range C^/R^ interaction even though the Hg(^Po) atom has no quadrupole moment. 
Using the conditions Ma + Mi, = M^ + M^ = Q, where Q is the total angular momentum 
projection and a good molecular quantum number, we get the matrix elements 

Yj{Sa,Si,,Ja,Jb,^) I -| \n D~5 



X 



{2-Ma + M'J.{2 + Ma-M'J. 

Ja 2 Ja \ I Jb 

-Ma Ma-M'a M'a ) \ Ma - 

24 

±(1 - 5a,b' 




{a\\Q\\a){b\\Q\\b) 



{2-Ma-M'a + n)!(2 + Ma + M'a-^)\ 

Ja 2 Jb \ I Jb 2 Ja 

-Ma Ma + M'a-n - MM I - -(M„ + M^-l]) M'a 



(15) 



where the dimensions of the matrix are given by max(f2 — J^, — Ja) < Ma < min{Q + J^, Ja) 
and max(fi — Jb, — Ja) < M'^ < mm{Q + Jb, Ja). 

The C5 coefficients are given by the eigenvalues of this matrix. The eigenvectors can be 
used to complete identification of the molecular state S = {Sa, Ja, Sb, Jb,^,'^,'R;i} where 
0<Q<4, 'n' = goTu, TZ = + OT—, and i distinguishes states of the same symmetry. In 
a few cases of double zero eigenvalues it was necessary to take the symmetric and antisym- 
metric linear combinations of the degenerate eigenvectors. 

The resulting C5 coefficients are given in Table [Tll It can be seen that of the total 
90 distinct molecular states, 37 have positive (repulsive) C5 coefficients, 32 have negative 
(attractive) C5 coefficients, and 21 have zero C5 coefficients. Of the 21 states with zero C5 
coefficients, 12 are "nonobligatory" in the sense that nonzero values would be allowed by 
the selection rules. 
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III. POLARIZABILITIES OF ATOMS AND MOLECULAR Cq COEFFICIENTS 



The dispersive van der Waals interaction between two atoms, Cq/R^ can be accurately 
calculated in terms of the dynamic polarizabilities of the two atoms. It can be approximately 
calculated in terms of the static polarizabilities of the two atoms, and a^- Only the latter 



are available for excited mercury atoms. The most frequently used approximation ^ 
the Cq coefficient in terms of the static polarizabilities is the Slater-Kirkwood formula 



for 



13| 



= ^ ^ -^""^ (16) 

2 ^/aa/ua + ^/ab/Ub 

where na and rib are the numbers of electrons in the outer shells of the two atoms. However, 
there exists relatively little experience with excited states. To validate the approximation 
and get some idea of its accuracy we sought similar systems where both accurate static 
polarizabilites and accurate van der Waals coefficients are available. Two apropos systems 
were found: (1) rare-gas dimers with both atoms in excited states and (2) zinc dimers 
with one atom excited (4^P or 4^P) and the other atom in the ground state. In both cases 
the Cg coefficients were calculated using the exact Casimir-Polder formula 1^ in terms 
of the frequency-dependent polarizabilities. The comparison of the Slater-Kirkwood and 



Casimir-Polder results in Table [TTTl suggests an accuracy o 



The ab initio values obtained by Rosenkrantz et al. 171] for the static polarizabilities of 



20%. 



the mercury atom in its 6s6p and states with S(Mi = 0) and n(Mi^ = 1) projections, 
neglecting spin-orbit coupling, are given in Table |Vl These values have been transformed to 
the J-Mj representation using the transformation matrix given in their paper, which was 
determined by fitting the experimental atomic term values [is!]. This transformation assumes 
that the radial functions of states with the same L, S but different J, Mj are essentially the 
same. No experimental values are available for the excited states, but their values for the 
ground-state of mercury, a = 36.2 calculated with a 2-electron relativistic effective core 
potential (ECP) and 34.3 calculated with a 12-electron relativistic ECP, can be compared 
with the experimental value, a = 33.91 ± 0.34 19|. The excited-state polarizabilities that 
we use for the triplet states were calculated with the 12-electron relativistic ECP, but the 
singlet excited states were calculated only with the 2-electron relativistic ECP. 

Rosenkrantz et al. 171] also did a similar calculation of the static polarizabilities of the 
excited zinc atom with a relativistic 2-electron ECP. We may get some additional idea of the 
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accuracy of the mercury polarizabilities by comparing these results with those of Elhngsen 



et al. 20|] for zinc, given in Table UTTl The results of Elhngsen et al. presumably should be of 
of higher accuracy since they account for core polarization; also, their value for the ground- 
state polarizability is in excellent agreement with the experimental value. This comparison 
is shown in Table IIVI The triplet and singlet excited-state polarizabilties of Rosenkrantz et 
al. are ~ 11% and ~ 78% larger, respectively, than those of Elhngsen et al.. This agreement 
is quite satisfactory for the triplet states, but less so for the singlet states. 

These comparisons suggest that the uncertainty in the Cg coefficients, which we present 
below, arises more from the static polarizabilities themselves than from use of the Slater- 
Kirkwood approximation. 

The eigenvectors x of Eq. f[T^ were used to calculate the polarizabilities for the atomic 
states forming each molecular state 

= x^(5)Ax(5) (17) 

where A is a diagonal matrix containing the polarizabilities of Table |V^b) on the diagonal. 
The corresponding Cq coefficient is then obtained using Eq. (fT6|) . These eigenvectors, which 
were obtained for the leading asymptotic interaction, are expected to still be a reasonable 
approximation when higher-order terms are included. The resulting Cq coefficients are given 
in Table IIIll They can be seen to vary significantly among the various states, with the 
interactions involving the singlet atom generally being the larger because of its greater 
polarizability. However, the comparison in Table IIVI suggests that the Rosenkrantz et al. 
polarizabilities for the singlet states may be significantly too large. If these polarizabilities 
are reduced by the same factor, 0.56, as for zinc, the coefficients would be reduced by 
a factor of about 0.66 for the molecular states arising from a ^Pi atom interacting with a 
'^Pj atom and a factor of about 0.42 for the molecular states arising from the interaction 
of two ^Pi atoms. At least for zinc, the uncertainty in the triplet atomic polarizabilities 
appears to be much smaller — the values of Elhngsen et al. being only 10% smaller than 
those of Rosenkrantz et al. Similar changes for mercury would reduce the Cg coefficients 
by about a factor of 0.86 for molecular states formed from two triplet atoms and by about 
a factor of 0.92 for molecular states formed from a singlet and a triplet atom. With both 
changes simultaneously, the effect on the Cq coefficients for a singlet and triplet pair are 
approximately multiplicative. 
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IV. ASSOCIATIVE IONIZATION CROSS SECTIONS 



The long-range potentials are very important, even dominant, for AI in low-energy col- 
lisions of two excited mercury atoms. For the present estimates of the cross sections, we 
adopt a model similar to that used in Ref. i.e., ionization is assumed to occur if the 
trajectory on the quantum-mechanical potential curve has a classical turning point smaller 
than both the black-sphere distance Rg and the crossing distance (if any) of the poten- 
tial curve with the Hg^ potential curve. The choice of -R^ = 4 A is discussed in Ref. 
The sensitivity to this choice will be examined below. The relevant data characterizing the 
short-range potentials, from Ref. is given in Table IVTl ^ For each potential curve, the AI 
cross section is taken to be 

a = 7r[min(6^, b^, 6o)]^- (18) 

The first two impact parameters depend only on the short-range potential. The largest 
impact parameter accessing the black-sphere distance is 

l-^V%. (19) 



E 

and the largest impact parameter accessing the crossing into the continuum is 

1-^)'%.. (20) 

assuming the collision energy E exceeds the potential energies V{R) at these points. The 
value of bo, determined by the long-range potential, is given by [see Appendix] 



where Ro is obtained by numerically solving 

Ro dV{Ro) 
2 dR 



+ ViRo) = E, (22) 



^ A clarification may be helpful. In Ref. [l| the cross sections ignoring the barriers in the ab initio potential 
curves (labeled "nb" there) were calculated neglecting the potential only at i? > 5 A. In the present 
calculations, the short-range potential is taken to be completely characterized by the values of Rs, T4, 
Rx, and T4 in Table IVTl (interactions not listed in this table have potential curves too repulsive for the AI 
reaction to occur). Also, note two typos in Ref. [ij: (1) In Table II, ^Pi + '^Pq has only one state of 0~ 
symmetry, and (2) in the caption of Fig. 3, 0~=dash dot. 
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with 



+ (23) 



for its largest real root. 

It can be seen in Table I VI I that most of the contributions to the AI cross sections at a 
collision energy of 0.00095 a.u. are limited by the long-range potentials, which determine the 
classical orbiting impact parameter. However, it should be kept in mind that the potential 
energy at the curve crossing (if any) and the black-sphere distance provide a necessary 
condition. For AI to occur, the collision energy must exceed the potential energy at these 
distances. Thus, without tunneling, which can be expected to be negligible except possibly in 
ultracold collisions, the cross section will vanish at sufficiently low energies if the interatomic 
potential energy is positive at either of these distances. At the collision energy of 0.00095 a.u., 
47 of the 90 potential curves do not contribute to AI for this reason. At still lower collision 
energies, the contributions that are energetically allowed are increasingly determined by 
the long-range potentials. At ^ 2 x 10~^ a.u., the non-zero contributions are entirely 
determined by the long-range potentials. 

We now examine the sensitivity of the cross sections to the parameters, in particular Ce 
and Rs, which are most uncertain. The classical orbiting cross section is simply obtained 
for pure power- law potential energy curves [see Appendix]: for V{R) = C^/R^ it is c'"orb{5) = 
¥ (^)^^^ and for V{R) = C^/R^ aoA{6) = Svr (^)^^^ With the dependence Cl^\ a 
reduction in Cq by a factor of 0.66 would decrease the cross section by only a factor of 
0.87, so the sensitivity is not too great. The uncertainty due to the parameters associated 
with the short-range potentials, in particular whether Vg and are positive or negative, is 
greater than this. 

For low-energy collisions, there is little uncertainty due to Rs, which is more-or-less 
arbitrarily chosen (see Ref. This does not necessarily mean that the magnitude of 

the autoionization width, which is presently unknown, is unimportant. It does mean that, 
accepting the black-sphere model, the value taken for Rg is not critical. If the width is large 
enough to saturate ionization in the energetically allowed region, its exact magnitude will 
be unimportant. If not saturated, the width will reduce the size of the cross section by a 
factor that approaches a constant at low energies. 

The total AI cross sections are given by statistically weighted sums over all the potential 
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curves arising from the pair of atomic reactants i = 2S+ip^ ^^^^ ^/ ^ 2S'+ip^^^ 

= (2J, + ir\2J, + ir\2 - 5,,0"' 5^(2 - SnM+I'^ (24) 

where the sum goes over the molecular states listed in Table [TTl for each reaction. This sum 
reflects the double degeneracy of the Q ^ states and the pairing of g-u states in the cases 
of nonidentical states. 

The resulting cross sections are shown in Fig. [H These cross sections are expected to 
be valid down to energies where the scattering is predominantly s-wave. The maximum 
contributing angular momentum is roughly 

L^,. = {2iiEY'%o. (25) 

where /i is the reduced mass and the impact parameter ho is given in the Appendix. For an 
attractive C^jP^ potential, L^ax ~ 847(-C5)^/^i?^/^°, and for a Cf^/R^ potential, L^ax ~ 
831(— Ce)^^^-^^^^ for Hg*+Hg* collisions. Thus the semiclassical description, which requires 
-Z^max ^ 1, should be valid above collision energies of ~ 5 x 10~^^ a.u. (temperature of ~ 2 
/iK) with typical C5 coefficients and above ~ 5 x 10^^^ a.u. (temperature of ~ 20 /xK) when 
all C5 coefficients vanish. 

Given the uncertainty in the short-range potentials, it is of interest to consider another 
approximation to the AI cross sections that might have validity at very low collision energies. 
Such an approximation can be obtained by considering only the long-range potentials, i.e. 

CTLR = Tihl (26) 

instead of Eq. [181 This approximation neglects any reduction due to repulsion in the short- 
range potential. These cross sections are shown as dashed curves in Fig. [H 

At very low energies the dependence of the cross section on energy tends to be a power 
law: E~'^/^ in cases with contributing attractive C^/B!' potentials and E~^/^ otherwise. At 
sufficiently low energies, potential curves with positive C5 coefficients will cease to con- 
tribute to the cross section, and potential curves with negative C5 coefficients will have this 
contribution dominate over any contribution from the Cq potential. The thermal average 

-<T) = a{E)Ee-'^l'^dE (27) 
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is particularly easy to obtain for a cross section having power law dependence. For 



a{E) = cE"", (28) 

a{T) = T{n + 2)a{kT). (29) 

where T is the Euler gamma function, with pertinent values r(— 1 + 2) = 0.9027 and r(— 1 + 
2) = 0.8935. 

The cross section for AI in the ^Pq + '^Pq collision is displayed dotted because it is a 
special case, subject to additional uncertainty. This cross section is presented as an upper 
bound, as are all the black-sphere cross sections, but, unlike the other cases, the true cross 
section is expected to be considerably smaller. In order for this reaction to occur at all in 
cold collisions, the Hg(^Po)+Hg(^Po) asymptote must be above the minimum in the Hg2' 
well. Though not all previous work has agreed on the sufficiency of the binding energy of 
Hg^, the most reliable experimental and theoretical results now agree that this well depth 
is greater than the required 1.10 eV (see Ref. [1] for discussion of this point). However, the 
single potential curve coming from the Hg('^Po)+Hg(^-Po) asymptote does not cross the Hg^ 
curve, so adiabatically AI still would not occur. The reason AI is nonetheless possible is 
that there is an avoided crossing of this curve with another curve of 0^ symmetry coming 
from the Hg(^Pi)+IIg(^Pi) asymptote. Their diabatic crossing probability is ~ 0.33, only 
weakly dependent on the collision energy as long as E is small compared with the potential 
energy of 0.013 a.u. at the crossing. An important point here is that the state corresponding 
to this potential curve has 93% quintet character and quintet states do not couple to the 
ionization continuum (only singlet and triplet states can be formed by coupling Hg^ to an 
electron, and the relevant coupling matrix element conserves spin). Thus the autoionization 
width is reduced by a factor of ~ 14. Even if the black-body model is valid for the other 
reactions, it is much less likely to be valid with an autoionization width this small. 



V. CONCLUSIONS 

The uncertainty of long-range barriers, which affected the results in Ref. [l|, has been 
resolved. The long-range 1/P^ interactions, which may be attractive or repulsive, have been 
accurately calculated. The 1/P^ interactions, which are generally attractive, have been 
determined less accurately, but still sufficient for useful estimates of the AI cross sections. 
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The only AI cross section for mercury that has been experimentally measured is for 
Hg(^Pi)+Hg('^Po)- The most reliable values are 99 and 64 A^jWith error bars of 25 A^, 
which come from two different analyses of the same experimen t |22l| (see footnote 5 of Ref. 

although there is an older published experimental value 23| of 460 A^. The present 
thermally averaged value for temperature 300 K is 83 A^, in agreement with the most recent 
experiment. This agreement tends to validate the black-sphere model, which, in principle, 

n 

provides an upper limit on the cross sections. Another experimental value [24] of 160 ± 40 



A^ has been published for the 



on experimental grounds 



22 



Hg 



25 



Po)+Hg(^Po) reaction, but has been generally repudiated 



261] and is inconsistent with the present upper limit on the 
cross section for this reaction. 

The cross section for AI in the Hg(^Po)+Hg('^Po) deserves attention because of the pro- 
posed use of ^Pq as an atomic clock state, even if the discussion is necessarily speculative. 
Compared with the other reaction pairs, the "strength" of this reaction may be reduced by 
three factors: a factor of ~ 0.33 due to the required curve hopping, a factor of ~ 0.07 due 
to the non- ionizing quintet component, and a factor due to the time spent in the continuum 
being fairly short. Only the first factor has been taken into account in the cross section 
displayed as a dotted line in Fig. [1] The other two factors could easily reduce the cross 
section by another order of magnitude or more. Better theoretical estimates will require 
determination of the relevant autoionization widths. 
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APPENDIX: CLASSICAL ORBITING CROSS SECTIONS 

The classical orbiting distance Ro and impact parameter ho at collision energy E are given 
by the solutions of the simultaneous equations, 

Kfr(i?) = E (A. la) 

and 
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where 

Ves{R) = V{R) + (A.2) 
for R and h. The corresponding cross section is 

CTorb = T^hl- (A. 3) 

For 

V(E^ = 

R^ 



V{R) = ^, (A.4) 



f^orb(n) = ^ TTd ' (^-5) 



n-2\ 2E 
so aorb(5) = T (^)^^^ and aorb(6) = Stt (^)^^^ For 

V(«)^|h-§. (A.6) 

the Eqs. (lA.ip can be solved numerically. The largest positive real root, assuming such 



exists, gives the desired result; otherwise, classical orbiting does not occur. 
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TABLE I: Quadrupole matrix elements of the excited mercury atom (Hg 6s6p) in atomic units. 



(a) Reduced matrix elements. 

^ ^ 

3Po -10.22 

8.218 -15.78 3.131 

-10.22 15.78 -15.38 2.860 

^Pi 3.131 -2.860 -29.67 

(b) Conventional quadrupole moments. 

^Po ^Pi ^P2 ^Pi 

3.00 -7.35 -10.83 



16 



TABLE II: Long-range interaction potentials for Hg(6s6p)+Hg(6s6p) with spin-orbit coupling 

(Hund's case c), in atomic units. 
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TABLE III: Comparison of Cq coefficients calculated by the exact Casimir-Polder (CP) formula 
using the dynamic polarizabilities with those approximated by the Slater-Kirkwood (SK) formula 
(Eq. fTUj) using only the static polarizabilities a. These results are presented to validate the SK 
formula for excited states. 



(a) Rare-gas dimers ( 


'^a = nb = 1 


in the SK formula). 
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315.8 


315.8 


4994 


4209 


Xe(6s, 3p2,2)+Xe(6s, 


'^'2,2), 4<7 


387.4 


387.4 


7138 


5719 


( n 1 7i n r* n itti pvt; (71 — 


nfe = 2 in the SK formula). 














Oib' 














CP'^ 


SK 


Zn(4s2, 15)+Zn(4s2, 


15), is 


39.12 


39.12 


282 


260 


Zn(4s2, 15)+Zn(4s2, 


3P), 3S 


39.12 


89.08 


435 


471 


Zn(4s2, 15)+Zn(4s2, 


3p), 3n 


39.12 


55.21 


370 


335 


Zn(4s2, i5)+Zn(4s2, 


ip), Is 


39.12 


396.74 


1139 


1258 


Zn(4s2, i5)+Zn(4s2, 


ip), In 


39.12 


116.51 


674 


567 



"From Ref. 
''Froin Ref. 
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H|. 

T^l; dynamic polarizabilities were adjusted to agree with Ref. '15l in the static limit. 
2^. 
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TABLE IV: Comparison of static polarizabilities. 



Polarizability 

Ellingsen et al. [20] Rosenkrantz et al. [17] Experimental [21] 



Zn(4s2, 
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TABLE V: Excited mercury atom (Hg 6,s6p) dipole polarizabilities in atomic units. 

(a) Ah initio calculation in L-S representation (from Ref. [17]). 

s n 

3p 128.38 58.77 

532.58 158.08 

(b) Transformed to J-Mj representation. 

Mj = Mj = l Mj = 2 

^Pq 81.97 

70.03 95.11 

105.18 93.58 58.77 

521.31 156.55 
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TABLE VI: Distances R (in A) and potential energies V (in a.u.) of the black sphere radius (subscript s), the crossing with the 

molecular-ion potential curve if any (subscript x), and classical-orbiting (subscript o). The corresponding impact parameters b 
(in A) are given for a collision energy 0.00095 a.u. Only the potential curves contributing to AI cross sections at this collision 
energy are listed. The cross section is determined by the smallest 6 (shown in bold face). 
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"With curve crossing tciken into account. 
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FIG. 1: Associative ionization cross sections for collisions of two mercury atoms in the indicated 
states. The solid curves take into account the long-range and short-range potentials (Eq. (jlSp ). 
The dashed curves take into account only the long-range potentials (Eq. (|26|) ). The cross section 
for Hg(^Po)+Hg(^-Po) is shown dotted as a special case (see text). 
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